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Abstract 
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operator realizations of [-D, A] = Id. 
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1. INTRODUCTION 



The idea of discrete physics and discrete space-time is a very old one. To the best of 
our knowledge, the oldest reference in a physical journal is [1] 1 and it has been reconsidered 
many times in the past (see [3] for an extensive bibliographical review on the subject) . The 
basic motivations for all those efforts were: (a) the presence of ultra-violet infinities in the 
standard quantum field theory; (b) understanding the origin of the fundamental length 
(mass) scale in the Einstein's general relativity theory. 

In recent years there has been a growing number of attempts to find an underlying 
discrete structure of space-time. At the same time interesting discrete structures have 
surprisingly emerged within originally continuous models. The most interesting examples 
can be found in string theory [4], Roger Penrose's spin network calculus [5], the loop 
representation of quantum gravity [6] (see also [7]), the Bekenstein black hole entropy 
problem [8, 9, 10, 11, 12] or the recent approach advocated by G. t' Hooft [13,14,15]. Also, 
it is interesting to mention in this context the old work by Aharonov et al [16], where the 
modular variables for the coordinate and momentum operators have been introduced to 
describe an infinite slit experiment with magnetic field. 

Increasing popularity of discrete models in theoretical physics has been stimulated by 
several successful attempts to discretize continuous models for technical reasons. This has 
advanced our understanding of discrete techniques. Here, the best examples are the lattice 
QCD [17] and the Regge calculus in classical gravity [18]. 

On the other hand, there have been several attempts to incorporate discrete space- 
time at the kinematical level and to investigate its physical consequences [19, 20]. This 
line of research can be also associated with the (/-deformations of the Poincare group [21]. 
Another interesting approach has been recently proposed by Mazur [22, 12], who has given 
additional physical arguments in favor of space-time discreteness. 

We assume a discrete coordinate space (see also [23]) and investigate its consequences 
for the Heisenberg commutation relations. This approach is purely kinematical, as the 
Heisenberg relations are, and does not depend on the details of the underlying dynamics. 
Such " unusual" realizations of the Heisenberg algebra have also been found in 2D gravity 
[24] . These representations are not unitarily equivalent to the Schrodinger representation. 
In particular, they cannot be exponentiated to the Weyl form of the canonical commutation 
relation. 

1 The first questions about the world geometry were already posed by Riemann in XlXth 
century [2]. 
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To motivate our discussion we first consider the standard quantum mechanical repre- 
sentation of the momentum operator: 

P , = -in ^ . (i.i) 

Hence, the coordinate discretization can be implemented through a derivative discretiza- 
tion: 

I - D " • (1 - 2 > 

where D Ax is a discretized derivative (difference) operator and Ax is the discretization 
parameter (usually a fixed coordinate spacing). Even though we will call D Ax a discretized 
derivative it should be noted at this point that initially this is just some difference operator 
acting, for example, on the space of smooth functions on R. However, such operators can 
be naturally restricted to act on the space of functions on Ax x Z. It is in this sense that 
we talk about a discrete derivative. 

We will consider a wide class of possible derivative discretization schemes of the form: 

+N 

D Ax (M + N+l) = D Ax = a k E L, (1-3) 

k=-M 

where au are real constant coefficients and M, N are integer indices corresponding to the 
lowest and highest non-zero terms (hence, ct_M ■ oln ^ 0). We take the coefficients ctfc to 
be constant as one would like to have the same definition of the discrete derivative at all 
points in the coordinate space. The expression will be called n-point for n=M+N+l. For 
convenience, let us define the shift operator E Ax as: 

El x f{x) = E n f{x) = f(x + nAx) . (1.4) 
One can easily check that this operator fulfills the following commutation relations: 

[E m , E n ] = , 
[E m , x] = m Ax E m , 

(1.5) 

[E m , x E n ] = m Ax E m+n , 

[E m , E n x] = m Ax E m+n . 

Several authors have used a purely imaginary Ax and ctfc. This case is not addressed in 
the present paper and our methods are not directly applicable to it. 
In addition, the formal continuum limit reads: 
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lim El x = 1. (1.6) 

Ax— >0 

In the simplest case: N = M = 1 the difference operator can be expressed by the Heine 
bracket [ ■ ] q defined for the dimensionless variables P, q as (see Appendix A for its basic 
properties) : 

( q p - q-P) 

H(P,q) = [P] q = ^_ q 9 _ 1} J , (1.7) 

where q = e 1 and P = —i-J^. This notation is widely used in the context of g-deformations 
and, in particular, to define the g-deformed Poincare group [21] (see also [22]). 

This paper is organized as follows. In Sec. 2 we investigate the basic properties of the 
difference operator (1.3). Sec. 3 is devoted to a construction of the Heisenberg algebra for 
the operator (1.3) and its conjugate operator X Ax when the latter is a bounded operator 
(finite series). In Sec. 4 we consider the case of unbounded operators. Sec. 5 introduces 
the momentum formulation of our problem and within this representation we complete our 
proofs in the unbounded case in Sec. 6. Sec. 7 contains the summary of our results. 



2. THE DISCRETE DERIVATIVE 

In this Section we shall investigate the general properties of the discretized derivative 
(1.3) and the corresponding momentum operator. We would like to be able to talk about 
the convergence of discretized operators to the ordinary ones when the parameter Ax tends 
to zero and whenever we can make sense of this limit we will call it the classical limit. In 
particular, we stipulate that for analytic functions: 

D Ax f{x) — ^f(x) (Ax - 0) , (2.1) 

holds pointwise in x. Consequently, by writing (2.1) in terms of the translation operator 
e kAx ik we obtain: 

+N 

D Ax f{x) = a * ekAx£ /(*) • ( 2 - 2 ) 

k=-M 

Applying the Taylor expansion to (1.3) at x = Xi we have: 
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+ N 1 d ff x .\ + N 

Dax f{xi) = f(xi) ^2 a k + Yy d 1 E a k k Ax + 

k=-M ■ X k=-M 

+N V- 6 > 

+ ... +_/(»)(x 4 ) V a k k n Ax n + ... 

Ul k^M 

From (2.3), we then get the following two conditions for the coefficients a k to fulfill (2.1): 

+N 

E «fc = o , 

(2.4) 

+N I 

E * «* = ^ • 

k=-M 

We may also sharpen the notion of a classical limit by specifying the accuracy of the 
discretization scheme. For example, the coefficients can be chosen to give the best fit to 
the ordinary derivative by requiring that all terms in (2.3)vanish up to the order Ax M+N . 
As a result, in addition to (2.4), the following set of algebraic equations for the coefficients 
must be satisfied: 

k=+N 

k n a k = for n = 2, 3, . . . , (M + N) , (2.5) 

k=-M 

and we get instead of (2.3) the following approximation for the derivative: 



D Ax f(x t ) = /'(*,) + { m + N+1)\ f (M+N+1) ^ X 



+M 

X 

k=-M 



a k k M+N+1 + o (Ax M+N+1 ) . 



(2.6) 



Eqs. (2.4), (2.5) comprise the set of M+N+l linear eqs. for M+N+l coefficients a k 
and they can be solved explicitly. The corresponding determinant is the Vandermonde 
determinant and its value is: det ||A|| = 1! ■ 2! • . . . • (2N)l ^ 0. Consequently, there exists 
a unique solution. Its form is given by the following: 

Theorem 1: The discrete derivative of the form (1.3) gives the best fit to the ordi- 
nary derivative, i.e. it has expansion (2.3) with the first non-vanishing coefficient for 
Ax M+N+1 , if and only if the coefficients a k are of the following form: 
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-1 

a 



n 1 m 1 
fc=i fc fc=i fc 



Ax 

(-l) fc+1 M! AH 



(2.7) 



Ax k (M + k)l (N — k)\ ' 
Higher order terms in the expansion (2.3) cannot vanish, as the system would become 
over determined. 

In the symmetric case, i.e. when M = N, Theorem 1 implies the symmetry between 
forward and backward terms. Indeed, in this case the coefficients (2.7) have the following 
simple properties: 

a k (Ax) = - a- k (Ax) = a- k (-Ax) , (2.8) 

Moreover, for a fixed k, 

lim a k (Ax) = (-l) k+1 -i- (MO). (2.9) 
JV^oo {/\x)k 

The meaning of (2.9) is explained in Appendix B. Moreover, eq. (2.8) implies the symmetry 
of the discrete derivative: 

D Ax = +D- Ax . (2.10) 

The operator D Ax satisfying (2.10) we shall call symmetric (in Ax). Later on we shall 
limit ourselves to this case as it preserves the parity symmetry. Discretization schemes and 
corresponding coefficients satisfying, in addition to (2.4), (2.5) and symmetric we shall call 
optimal. In fact, in the special symmetric case N = 1 (2.8) implies that the derivative is 
optimal and we obtain the well known central difference scheme. The numerical values of 
the optimal a-co efficients for the lowest (2A^ + l)-point schemes are given in the Table 1 
below. 

In analogy to (2.2), the optimal derivative can be rewritten as: 

D Ax f(x) = 2 J2 a k smhikAx—j f(x) . (2.11) 
k=i ^ x ' 

In particular, for the optimal 3-point (= central difference) discretization scheme we have: 



B - /(*> = ^ -h m = rib sin A *l ) ^ ■ < 2 - 12 > 
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Table 1. Values of the lowest optimal a-coefficients. 



N = M 




«2 


«3 


Ct4 


«5 «6 


1 


1/2 










2 


2/3 


-1/12 








3 


3/4 


-3/20 


1/60 






4 


4/5 


-1/5 


4/105 


-1/280 




5 


5/6 


-5/21 


5/84 


-5/504 


1/1260 


6 


6/7 


-15/56 


5/63 


-1/56 


1/385 -1/5544 



The representation (2.12) is exactly the derivative represented by the Heine symbol (1.7) 
([21, 22]). 

Originally, Heisenberg introduced his matrix mechanics by postulating that mathe- 
matical operations in classical equations of motion should be reinterpreted (that has led 
to the correspondence principle). In this spirit Dirac has postulated that the "quantum 
differentiation" must satisfy the additivity and Leibniz rule. From this he has obtained 
the fundamental (Heisenberg) commutation relations for coordinates and momenta [25]. 
We would like to proceed in an analogous way. However, in our case not all general rules 
for derivatives are satisfied by (1.3), subject to (2.4). The basic properties of the ordinary 
derivative are listed below along with some comments on how those properties change if 
one uses D Ax instead of 4-: 

1. Additivity. Ordinary and discrete derivative are both linear : 

Dax [f{x) + g(x)] = D Ax f(x) + D Ax g(x) . 

2. Leibniz rule. The rule for differentiation of the product is not satisfied exactly by (1.3). 
The deviation is of order O (Ax) and for the optimal scheme of order O (Ax 2N ). 

3. Derivative of the composite function. The same result as for the Leibniz rule holds for 
the derivative of a composite function: D Ax [f o g(x)] = D Ax f(g(x)) x D Ax g(x) + 
0(Ax 2N ). 

4. Derivative of monomials. Discrete derivative of the monomial x n is identical to the 
ordinary derivative if the discretization scheme is optimal and the exponent n < 2N. 
In general, assuming (2.4) and (2.8) we have: 
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n 



Dax(N) x n = n x 



n-l 




Ax 1 X 



(2.13) 



5. Heisenberg commutation relation. This will be the subject of the following sections. 
We shall present a construction of the discrete counterpart of the Schrodinger couple x 
and ^ satisfying a discrete analog of the continuous Heisenberg commutation relation: 
[-^,x]f(x) = f(x). Naively, that is, by replacing only with D Ax and using (1.3), 
we get: 



For the optimal discretization we will construct a position operator called X which 
together with D/± x forms a conjugate pair, that is the relation [Dax, X] = I holds on 
a dense domain in a Hilbert space. 
6. Hermicity. In the symmetric case we have a natural Hermitian conjugation induced 
from the relation (E^ x )^ = E^™ The operator D/± x , like the operator ^ is formally 
anti-Hermitian if M = N and ak(Ax) = —a-k(Ax) and we have: 



A more detailed analysis will be given in the subsequent sections. 

3. LIE ALGEBRAIC DISCRETIZATION 

The notion of the Lie algebraic discretization has been introduced in [26]. It can 
be described as a special case of representation theory of the enveloping algebra of the 
Heisenberg algebra in which at least one of the generators is a difference operator. As a 
result, some linear difference equations have a representation theoretic meaning similar to 
the D-module interpretation of linear differential equations [27], in particular in Sec. 6. 
The most important application seems to be to describe quasi-polynomial solutions to a 
certain class of difference equations. On the other hand the whole procedure has certain 
appeal of deja vue in the context of basic quantum mechanics and from this point of view 
one can regard the authors' concept of the Lie algebraic discretization as a quantization 
procedure. 

Following [26], we start with (1.1) and we look for a conjugate operator X = X Ax 
that satisfies: 




(2.14) 



k 




(2.15) 



8 



[ D Ax , X ] = 1 , (3.1) 
with the additional condition ('classical limit'): 



lim Iat = x . (3.2) 

Ax^O 

In [26], only the special cases of the forward and backward Euler schemes that are not 
symmetric under reflections (2.10) have been discussed. This can be improved in two ways. 
Firstly, we can start with the more general discretization scheme (1.3). Secondly, we can 
require that -Dax be formally skew-Hermitian whereas X Ax be formally Hermitian relative 
to some * anti-involution to ensure a natural physical interpretation. 
In [26], it has been shown that for derivatives in the form: 

DiJ(x) s f±±^±zlM , D -J (x) s M , (3 .3) 

we get the Heisenberg algebra for D + and x(l — AxD~): 

[D+, x(l — AxD~)] = 1 . (3.4) 

However, as has been mentioned in Sec. 2, we are not satisfied with the Euler forms of 
the discrete derivative (3.3), as these schemes are not optimal and symmetric with respect 
to reflections 2 . In particular, this form of the discrete derivative does not seem to be skew- 
Hermitian with respect to any natural Hermitian structure. Also, the central difference 
scheme (and optimal, in general) was shown to enable a consistent formulation of the 
action principle and the Euler-Lagrange equations in the classical theory [28]. Another 
possibility would be that the pair D Axi X is a pair of conjugate operators (annihilation 
and creation operators). But here again, no natural Hermitian structure on the Hilbert 
space of states seems to induce that. 

First, we take up a question of generalizing the results of [26] to a wider class of 
discrete operators. As the starting point we generalize the form of the coordinate operator 
used in [26], stipulating the following Ansatz : 

M 

X EE & X ^ > ( 3 - 5 ) 

k=-N 

2 It can be shown that for classical dynamical systems the optimal schemes are closer to 
continuous derivative with respect to the stability analysis and they preserve the Hamilto- 
nian structure of the phase space, as well [A.Z. Gorski and J. Szmigielski, in preparation]. 
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subject to the condition (continuous limit): 



M 

E fa = 1 • ( 3 - 6 ) 

k=-N 

In fact, one can add to (3.5) two additional terms without changing the conclusions. Thus 
X might be taken to be of the form: 

M M M 

X = fa x E " + E fa E " x + E fa E k . (3.7) 

k=-N k=-N k=-N 

The last term commutes with the derivative D& x and does not change commutation rela- 
tion (3.1). It changes however condition (3.6) replacing it with: 

M M 

]T (p k + p' k ) = i, E fa = °- ( 3 - 8 ) 

k=-N k=-N 

The second term in the right hand side of (3.7), due to (1.5), has the same commutation 
relation as (3.5) leading to the substitution: 

Pk — Pk+P' k ■ (3.9) 

From now on, we restrict our attention to (3.5). From (3.1) we get the following conditions 
for the P l s: 

M 

E fa = 1 > 

k=-N 

n (3-10) 

Ax ^ a k (3 n - k k = 5 n0 , 

fe = -M 
-N<n-k<M 

where n goes from -(M + N) to +(M + N). Here, (3.10) is a set of (2M + 2A^ + 2) eqs. for 
(2M + 2A^ + 2) coefficients a k , f3 k plus two additional equations (2.4). Hence, the system is 
formally over deter mined. Let us start with the second equation above (set of 2iV + 2M + 1 
eqs.) and rewrite it as two separate sets, each of size M + N + 1, corresponding to n < 
and < n respectively. For simplicity we introduce matrix notation in the (M + N + 1)- 
dimensional space of coefficients. Let us define vectors: 
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p = 



P-N 
P-N+l 



b = 



- Pm 

and two (N + M + 1) x (N + M + 1) matrices: 



1/Ax 




A_ = 



N(*n ■ ■■ ••• -Ma_M 

(N-l)a N -! ... -Mol-m 



■Ma. 



■M 











A + = 



Ncxn —Ma-M 

A^ctiv ••• (-M + l)a_ M +i 



(3.H) 



(3.12) 



L ... Na N 

Now, our equations can be rewritten as two sets of equations for the coefficients Pk in the 
following simple form: 



A_ (3 



A+P = 6 , 



(3.13) 



where the first equations in both sets are identical. The determinants of the systems (3.13) 
are: 



detA_ = (-M «_ M ) M+N+1 , detA+ = (N M+N+1 a N ) M+N+1 . (3.14) 

and (3.13) can be solved explicitly. Indeed, there are three possibilities to consider: 1. 
both determinants are nonzero (i.e. M ■ N ^ 0), 2. one of the determinants is zero, 
then either M = or iV = 0, 3. both determinants are zero (M = N = 0). In the 
last case the derivative expansion (1.3) consist of one nonzero term (cto-Eo) only and the 
second condition of (2.4) cannot be fulfilled. In the first case, using Cramer's rule, we 
get two different solutions to (3.13), hence the system is inconsistent. The only nontrivial 
possibility is to have the second case. Let us assume that ^ 0. Then from the second 
equation of (3.13) we get the following (unique) solution for /3's: 

P-n = AT 1 A , P-N+l = ... =Pm = 0, (3.15) 
Naj^ Ax 

and the first equation of the set (3.10) implies: /3-n = 1 and «at = 1/NAx. Now, 
consistency of the first eq. of (3.13) implies that the only non-zero ck's can be: «o and 
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ql jv • Finally, eqs. (2.4) are satisfied as well if: a = —a^. Identical solution can be found 
for a-M 7^ 0. Hence, we have got the following 

Theorem 2: For the discretized derivative operator (1.3) satisfying (2.4) and the coor- 
dinate operator in the form (3.5) with condition (3.2) to satisfy Heisenberg commutation 
relation (3.1) the coefficients ctk^ftk must be either in the forward Euler form: 

Qfc = a k = ~ N \ x Sk0 + N \ x SN ' k ' Pk = Pk = +o~-N,k , (3.16) 
or in the backward Euler form: 

1 1 

ak = a k = ~M^ 6k0 - M~A^ 6 ~ M ' k ' Pk^PL 3 = +$-M,k ■ (3.17) 



Thus, for X of the form (3.5) the only way to satisfy (3.1) is to take Da x as the Euler 
derivatives (3.3). In this case the coordinate operators are: 

X B ' F = x (1 - Ax D+ A J = x Eg . (3.18) 

This proves that the Turbiner-Smirnov discretization scheme is unique under conditions 
(3.1), (3.2), (3.5). Hence, an optimal derivative with the discretized coordinate operator 
of the form (3.5) cannot satisfy the Heisenberg commutation relation. In fact, this is 
impossible for any symmetric discretization scheme with finite series in (3.5). 



4. TRANSCENDENTAL OPERATORS X 



In the previous Section we have limited our considerations to the coordinate operators 
(3.5) in the form of a finite series. However, there is still a possibility of having the 
Heisenberg algebra with the discrete derivative of the form (1.3) and an infinite series 
for X. We shall investigate the algebraic aspects of this problem in this section and a 
complete, rigorous solution of the problem will be given in Sec. 6. We should mention that 
the operator X is always unbounded and eventually a special care needs to be exercised 
when dealing with any expressions involving X. In this section we study only formal 
aspects of our problem leaving aside important questions like the domain of the definition 
of X or in what sense we understand expressions below involving infinite sums. We will 
come back to those questions in Sec. 5 and 6. 

In the present section we postulate that our coordinate operator has the following 
form: 
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X — x (3 m E% x , 



(4.1) 



k= — oo 



and we have an infinite set of linear algebraic eqs. (3.10) for n = 0, ±1, ±2, ±3, ... . In 
addition, the continuous limit condition (3.6) is tentatively taken to mean: 



+oo 



E 0* = 1 • 



(4.2) 



fc= — oo 



Eventually, we will have to regularize this condition since in many examples a direct 
evaluation of this limit does not make sense. 

Defining infinite dimensional matrices and vectors analogous to (3.11), (3.12): 



P+ 



P-N 
P-N+l 
fl-N+2 



P- 



Pm 
Pm-i 

Pm-2 



1/Ax' 





(4.3) 



A+ = 



Na N ... -Ma-M 
Net jv ... -Ma-M 



-Ma-M ••• Na N 

-Ma-M ••• Na N 



(4.4) 



formally, we get, instead of (3.13) the following two sets of eqs.: 



A + (3+ = b , A- (3- = b . (4.5) 

In both sets the first equation is identical and relates at most (M + N + 1) coefficients 
(3 m . Each next equation involves one more (higher order) coefficient that can be explicitly 
computed. Hence, our system has (M + N) free parameters. In the special case, when 
M ■ N = the first or last element in each row in (4.4) vanishes and we have (M + N — 1) 
free parameters. This implies for the kernel of operators A+, A-\ 



D = dim ker A + = dim ker A- = M + N - 5 M n,o ■ (4.6) 
Hence, we have proven the following 

Theorem 3: For operator X in the form of infinite series (4-1) there is (M + N — 5mn,o)~ 
parameter family of formal solutions to the Heisenberg relation (3.1). 
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Remark: We do not claim at this point that any of the formal solutions satisfy (4.2). 

In the case M ■ N ^ eq. (4.5) imply one algebraic condition for the coefficients 
{/3-n, ■ ■ -,/3m}- 

+N 

£ k ak ^ = Ax~ ' (47) 

k=-M 

and the higher coefficients fa (k > M and k < —N) can be uniquely determined from: 

-1 +N 
^ = MoZ^ S k a « » 

fc - M+1 (4.8) 



P-(N+i) = T7^- kak P-(i+k) > 



JV fc=-M 

where i = 1,2,... with additional condition (4.2). 

In the remainder of this section we shall assume that the coefficients a k are optimal 
and the coefficients fa have the symmetry: 

fa = fa k ,-N <k<N. (4.9) 

This implies the symmetry of the discrete coordinate operator with respect to the Ax 
reflections: 

X Ax = X. Ax , (4.10) 

in analogy to the discrete derivative operator (2.10). 

Indeed, now we have N + 1 /3-parameters {fa, fa, fa, ■ ■ ■ , 0n} subject to : 

N 

k=i 

The higher order coefficients {fa for \k\ > N) are uniquely determined from: 



_1 N i Ar " 1 
(3 N+i = — y^ k a k fa-k - -t-j y^ka k (3 i+k , z = 1, 2, . . . . (4.12) 

Using induction on i and (4.7) we get the claim: 



P-(N+i) = ftN+i 
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(4.13) 



We remark that the normalization condition (4.2) can be rewritten as: 



+00 

2 = (4.14) 

k=i 

i.e. we are left with (N — 1) free parameters (3. All considerations in this section were 
rather formal as we deal with infinite series whose convergence is unclear. In fact, as we 
shall see below the series might diverge in the usual sense. We discuss in details a concrete 
example to illustrate this point. 

In the case N = 1 (4.12) implies the following condition for the /^-coefficients: 

A+2 = Pi , i = 0,1,2,... . (4.15) 

Hence, we have two free parameters: Po and (3±. The latter can be determined from (4.11), 
giving: 

Pi = +1 • (4.16) 
This provides an exact formula for all odd coefficients in the following compact form: 

p2k+i = (-1)" , fc = 0,l,2,... , 

(4.17) 

P-(2fc+l) — P2fc+1 7 

and for even coefficients we have: 

p 2k = (-l) k Po ■ (4.18) 

It is clear from (4.17) that our series is indeed divergent and we must redefine the classical 
limit to make sense of the summation in (4.1). There is no unique way of doing that. We 
will say more about this in the next section. 

At the moment, we are still left with one free parameter, Po, and one equation, the 
normalization condition (4.14). As a special case we can put 

Po = . (4.19) 

Hence, a formal solution for the discretized coordinate operator, X, in the N = 1 case is 
(4.1) with P~co efficients given by: 

fok = 0, VA; 
P 2k+1 = (-l) k , k = 0,1,2,... , (4.20) 
P-(2k+i) = P2k+i , k = 0, 1, 2, . . . . 
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Now, this solution can be compared with the solution found by Frappat and Sciarino 
[29]. Their solution, after restricting it to 1 dimension, reads: 



1 1 



(4.21) 



X = x 



E + E- 1 E + E- 1 



x . 



There is a formal analogy between their solution and (4.20). To see that analogy one 
expands the first term in (4.21) at E = and the second one at E~ x = 0. There is, 
however, an important difference between the approach of the present paper and that 
of [29]. In their approach, the discretization parameter Ax is purely imaginary. As a 
consequence, the shift operator E is formally Hermitian and positive definite, while the 
momentum operator corresponding to our D& x is unbounded. Moreover, since the shift 
is in the imaginary direction it is hard to see what this has to do with a discrete space- 
time. By contrast, we are treating the shift as a unitary operator, thus retaining the real 
discretization parameter, which, in turn, ensures that the momentum is bounded. 

5. MOMENTUM FORMULATION 

In this section we reformulate some of the results obtained earlier as well as clarify 
some of the points left obscure in our discussion of equation (4.21). 

The basic idea is to realize the algebra of shift operators as an algebra of multiplication 
operators on the Hilbert space of square integrable functions on the circle. One can also 
view this realization as the momentum space counterpart of the content of the preceding 
sections. Throughout this section Ax = 1. 

At the core of the idea of discretization is that in all operations one proceeds by 
performing a sequence of discrete steps of size, in this section, 1. This calls for bringing 
the Fourier transform into the picture. To this end we will first reinterpret (1.5). Our 
setup is as follows. We consider the unit circle: 



We use the following parametrization of the circle S 1 : z = exp z6>, 9 G M (mod 2ir). 
Thus a function / e L 2 (5' 1 ) can be written as: 



S 1 = {zeC 

and the Hilbert space H = L 2 (5' 1 ). We recall 



z\ = l} 



(5.1) 




(5.2) 



16 



f(e^) = a n e %nd , (5-3) 

and a n are the Fourier coefficients of /. Since the shift operator E defined in (1.4) cor- 
responds to the elementary step of unit size we postulate that E = T 2 , where T z is the 
operator of multiplication by z. We will say more about multiplication operators towards 
the end of this section. Till further notice we also suppress writing T. Thus we set 

E n = z n , x = -z 4- , neZ 
dz 

It is routine to check that (1.5) hold. Furthermore we define 

N 

D(z) := a * zk ■ ( 5 - 5 ) 

k=-M 

In other words the difference operator becomes now a multiplication by a Laurent 
polynomial in z. We might observe that the operator x above is one of the generators 
of the Virasoro algebra (with the trivial central term), called, L . Now we reinterpret 
equations (2.4) and (2.5). We see that (2.4) goes into 

D(z = 1) = 0, (L D)(z= 1) = -1 . (5.6) 

Similarly, (2.5) goes into 

{L%D){z= 1) = , for n = 2,3,...,(M + iV) . (5.7) 

Now, it is immediate that the coordinate operator X should be sought (see (3.4)) in the 
form 

x = ~ z £ «■ • 

n= — oo 

where 

oo 

Pn = l- (5.9) 

n= — oo 

Let us now, at least formally, introduce the function h(z) = J2°^oo^nZ n . From the above 
condition we get that h(l) = 1. Moreover X can be written 
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(5.4) 



(5.8) 



X(z) = L h(z) = -zh(z)4- ~ zh'(z) . 

dz 



(5.10) 



As was already explained earlier, we can in fact ignore the second part and look for X in 
the form: 



X(z) = -f(z) 



d_ 

dz 



/(I) 



(5.11) 



Thus X is formally an element of the Virasoro algebra. We want that 



[D(z),X(z)\ 



1 



(5.12) 



where D(z) is a function (Laurent polynomial) defined previously. We therefore obtain 
that: 



f(z) D'(z) = 1 , 



(5.13) 



which formally implies 



m = 



D\z) • 



(5.14) 



We note that the condition f(l) = 1 is automatically satisfied as a consequence of the 
second condition imposed on D. However, it is important to underscore the fact that the 
latter expression can only be understood formally at this stage of our analysis. Yet, there 
are some special cases when we do not expect to have any difficulties. For example, if 
D(z) = A + Bz, that is M = 0, N = 1. Then we impose the following conditions: 



£>(1) = 0, £>'(1) = 1, (5.15) 

which gives D{z) = z — 1 and, consequently, f(z) = 1. Thus we get the pair D(z) = 
z — 1, X(z) = — which clearly satisfies (5.12). Using (5.10) we get that in this case 
h(z) = \. Hence = 1 and ft = 0,i ^ —1. Going back to the difference operator 
realization we get that D = E — /, X = xE~ x which is the pair found in [26] for Ax = 1. 
We can contrast this example with that of the symmetric discretization scheme for which 
D(z) = \{z — z~ x ). In that case f(z) = which is not bounded on the unit circle 

and consequently f(z) is not a multiplier (that is the multiplication by / is not a bounded 
operator on L 2 (5' 1 )). This example is explained in detail in the next section. 
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However, there are many cases for which / is a multiplication operator on the circle. 
We present now a general construction of X for which / is a multiplier. One can view this 
construction as a generalization of [26] . 

First, we review a few elementary facts from the theory of doubly infinite Toeplitz 
operators. Initially, we can consider a vector space of sequences a = {a n } ne z such that 
Ylmei \ a n\ 2 < 00 • This space, called £ 2 , is isomorphic to L 2 (5' 1 ), the fact well known from 
the theory of Fourier series. On £ 2 we introduce a linear operator T c : 

a m ~~ ^ ^ ] c m-n a n ? (5.16) 

where c = {c n } nez is a sequence of numbers. It is known, that T c is a bounded operator 
on I 2 iff c is a sequence of the Fourier coefficients of an essentially bounded function on 
S 1 . We recall that (f)(9) is an essentially bounded function on S 1 if it is bounded almost 
everywhere there. The isomorphism between I 2 and L 2 allows one to study the Toeplitz 
operator T c as acting on L 2 (5' 1 ). We will say that (f)(6) is a symbol of T c if under this 
isomorphism T c acts on L 2 (S l ) as: 

(T c f)(9) = (f)(9) f(9) , feL 2 (S 1 ) . (5.17) 

Because of this relation we often write to denote the Toeplitz operator with symbol 
(f). We will need in this section one more result from the theory of Toeplitz operators. 
The result below gives a complete description, in terms of symbols, of invertible Toeplitz 
operators. 

Lemma 1: is invertible if and only ifl/(f> is essentially bounded. If this holds, T^ 1 = 
T 1/4> . 

Now we apply this lemma to our problem. Let us first consider a Laurent polynomial 
g(z) such that Res z= og = and g is nowhere zero on S 1 . We can normalize such a g to 
satisfy g(l) = 1. Now we define 

r e 

D(e l °) := i / e lt g(e lt ) dt , (5.18) 
Jo 

where we used dz = ie l0 d9. Observe that since there is no z~ x in g, D(z) is again a 
Laurent polynomial satisfying (5.6). Now we set f(z) = j^j. We see that f(z) is bounded 
on S 1 so, by the lemma above, T/ is then the bounded inverse of T g . Now we define the 
operator X(z) = —Tf-i^: C ca (S 1 ) — > C°°(5' 1 ). By exactly the same computation as the 
one leading up to (5.5) we obtain that on C°° C L 2 (S' 1 ), X(z) and D(z) so defined satisfy 
the Heisenberg commutation relation as well as (5.6). 
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Our original setup has been formulated in terms of the £ 2 space. Now we would like 
to return to it. The map which maps back L 2 (5' 1 ) to £ 2 is given by the Fourier series 
method: / — > {f(n)} nez , where f(n) is the n-th Fourier coefficient of /. We recall that 
using (5.10) we can describe the coefficients (3 n as: 

p n = f{n+l), neZ. (5.19) 

Example: g(z) = a + bz~ 2 , \a/b\ 7^ 1. The condition g(l) = 1 gives a + b = 1. Thus 
D(z) = az + (1 — 2a) + (a — |a/(l — a)\ 7^ 1. We have two cases to consider 

depending on whether |a/(l — a) | < 1 (Case 1) or \a/(l — a)\ > 1 (Case 2). 

Case 1: f{z)=z 2 az , + \ 1 _ a) . 

We observe that / is analytic inside the unit circle 5" 1 and thus its Fourier expansion 
coincides with its Taylor expansion around z = 0. We get 

/ >) = ( ( " 1) i-: 2)/2 (l^) (n " 2)/2 ^ ifn>2andniseven; (5 _ 2Q) 

1 0, otherwise. 

Thus 

Pn = !^^(^-J n - 1)/2 , if n > 1 and n is odd; (5 . 21) 

1 0, otherwise. 

Case 2: f(z) = 0+(1 _ 1 a)z - 2 , |a/(l - a)| > 1. 

We observe that / is analytic outside of the unit circle 5" 1 and thus its Fourier expansion 
coincides with its Taylor expansion around z = 00. We get: 

f{n) = { I ^ ! ( i ^) ( " n)/2 ' if n < and n is even; (5 22) 
1 0, otherwise. 

Consequently, 

P n ={ i - a v—( i ?) <_ " _1)/2 ' ifnS-landnUodd; (5 . 23) 
1 0, otherwise. 

This example illustrates very well the fact that the optimal discretization is special 
even from the point of view of the Heisenberg commutation relations. Indeed the case 
excluded from the example is a = 1/2 which is precisely the optimal case. We can get, 
however, some information about this case making a approach 1/2. This allows us to 
discuss the formal solution (4.20). One way of interpreting this solution relies on a regu- 
larization of D which essentially amounts to moving the zeros of D' off the unit circle. We 
now present a simple example of such a regularization. First we set a = \ — e, e > 0. We 
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observe that | | < 1, so, we are dealing with the first case above. Directly from (5.21) 
we get: 

/j n (e) = ( ( "S-r (^fe) ( "- 1)/2 , ifn>landnisodd; {5U) 
[ 0, otherwise. 

It is interesting to note that 

a fn \ /2(-l)( n " 1 )/ 2 , if n > 1 and n is odd; oc .^ 
^ (0) = \0, otherwise. (5 ' 25) 

This is a formal solution to (4.7) and (4.8), with no symmetry condition imposed on /3's. 
It has exactly the same status as the symmetric solution (4.20). Thus we see that the 
formal solutions of Section 5 can be justified through an appropriate regularization. In 
addition we have to properly interpret (4.2) (classical limit). It makes sense to consider 
this expression for e ^ 0. Then, as it is easy to check, 

oo 

Pk(e) = l. (5.26) 

k= — oo 

The same expression does not make sense for e = 0, clearly indicating that one cannot 
interchange lim e ^o with 

The final item we would like to discuss in this section is the question of Hermicity or 
rather lack thereof. In our formulation the pair D and X may not necessarily consist of 
(even formally) Hermitian or adjoint elements. This is plain, for example, if one looks at 
the pair D = E — 1 and X = xE~ x or even more general pairs found above. One could, 
however, make this pair adjoint to each other, thus turning D and X into a pair consisting of 
a creation and annihilation operators, by constructing a proper representation space. One 
possibility would be to consider a Fock representation, which carry an, essentially, unique 
Hermitian form. Thus, for the last example, the completion of the space span{e n = 
X n l,n > 0} with respect to the inner product < e m ,e n >= S m ^ n n\ gives us a Fock 
space on which D acts as an annihilation operator and X acts as a creation operator 
respectively, and both act as difference operators. Another possibility is to consider a 
coherent representation, i.e. such which is generated from the vector v : Dv = Xv, A G C, 
whose special case is the Fock representation obtained for A = 0. In the case of D = E — 1 
the vacuum state v is a quasiperiodic function satisfying v(x+l) = (X+l)v(x). One obtains 
the representation space as the span{e n = X n v,n > 0} and the unique Hermitian form is 
defined by< v,v >= 1,< e m ,e n >=< u, D m X n v > where in the last expression one moves 
X to the left and D to the right making use of (5.12), the formula < it, X n w >=< D n u, w > 
and Dv = Xv. 
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In the next section we study the optimal discretization for which we show that one 
can construct pairs of Hermitian (self-adjoint) operators D and X. 

6. MOMENTUM FORMULATION APPROACH TO OPTIMAL 
DISCRETIZATION 

Now we take on the case of the optimal discretization. Thus M = N throughout this 
section. Also, in this section we are interested in self-adjoint pairs X, P relative to the 
standard inner product on L 2 (S 1 ). For convenience we multiply D from previous sections 
by 1/i and denote the resulting function by D^. 

The problem of classifying the pairs of operators satisfying Heisenberg commutation 
relations under an additional assumption that one of the operators be bounded was con- 
sidered in [30]. Since our difference operators are bounded, we are studying a special case 
of that classification, included in [30], in particular in Theorem 8.5 therein. The approach 
below is direct, however, and can be used as an introduction to a more encompassing 
treatment of [30]. Moreover, we have a different physical motivation, our interest lies in 
concrete bounded operators explicitly given by difference operators, reflecting the under- 
lying assumption that the coordinate space is discrete. We recall some basic definitions 
from [30] . Let H denote a separable Hilbert space over C. 

Definition: Let P, Q be operators in H, and O a dense subspace of H. We call (P,Q) a 
conjugate pair on Q iff 
(Kl) P and Q are symmetric, O C D(P) n D(Q), 

(K2) Pfi C O, Qfi C O, 

(K3) P = Pjtt, Q = Q/n, 

(K4) Q is bounded. 

In this section we essentially show that our X and form a conjugate pair in the above 
sense. 

We start with 

Lemma 2: Let Dn(z) be optimal, then D' N (z) has two simple roots on the circle S l . 

Proof: Since z = e ld , < 9 < 2n, we want to show that j§D N = D' N {6) has exactly two 
zeros for 9 e [0, 2ir). 

To derive D' N {9) we first compute Dn(9). We get that 
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N 



Hence, 



^-r — ttt TTTrrr re, sinkd • 

fc (iV + /c)!(iV- &)! 



N r 7vn 2 



or, in terms of z: 

D *W ~ k ^ N (N + k)\(N-k)\ 

fc^O 

_ (-l)^(iV!) 2 „ " , ^ k (2N 



z = 



(2N)\ ' ^ > \ N 

fe^O 



k = -N V 



z k = 



(_l)^+i(JV!) 



2 



z- N 



(1 - ,r-(-i)"( 2 ^)^ 



(2AT)! 

Thus it suffices to show that the polynomial (f — z) 2N — (— 1)^ (^jj^j zN has exactly two 
zeros for z <E S 1 . Indeed, we obtain: 



'This equation lias two solutions, provided: Z \l [ ^jy ) < ' ■ ^ — ' ■ 



Indeed, for iV = 1 we get |-2 < 1 and there are two roots, 9 = ^ , §7r. Now, we proceed 

t.n- I 



2iV 

by induction. First we observe that the above inequality is equivalent to: I AT I < 4^. 



We already proved it for iV = 1. We note that: 

2iV + 2\ /2iV\ (2iV + l)(2iV + 2) / 2N \ 2(2iV + 1) 



N+l J \N J (N + l) 2 \N J (N + l) ' 

and we use the induction hypothesis to prove the claim. 

To prove that the zeros are simple we show that ^[(1 — z) 2N — (— 1) N ^^y^ 

not zero there. Using that (1 — z) 2N — {—1) N ( "^T ] z N = we obtain that this derivative 



z N ] is 



x N , 

can be zero if and only if 2z 2 — 2z — 1 = 0. Since this equation has two real zeros not equal 
to ±1 we conclude that the zeros of D' N (9) are simple. ■ 

Now, we turn to the problem of the localization of those zeros. We will show that 
both zeros move towards 9 = n as N — > oo. First we prove three technical lemmas. 
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Lemma 3: 



(1+2*) < 11 



2(77+1) 



N > 1 . 



N+l 

_J_ ) _ 2N 2 +5N+2 j /-, , iwi , 1 \ _ 2jV 2 +5iV+5 
AT+1^ 2AT(AT+1) dlm l ir jyA 1 ^ 2 (AT+i)J 2JV(iV+l) 



Proof: Consider (1 + ^)(1 + 
Hence, (1 + ^)(1 + 

jv+t) <(! + *)(!+ 2 (at+i) ) which implies the claim. 
The next lemma improves on the estimate used in the course of the proof of Lemma 

2. 

Lemma 4: 



2N 
N 



< 4 



N 



N > 1 



(1 + *) 

Proof: The proof goes by induction on N. The only nontrivial step is to observe that: 
2(N+1)\ (2N\ 4(1+^) 
N + l J \ N 



which, by the induction hypotheses, implies: 



N+l 



2{N + 1) 
AT+ 1 

Now, the claim follows from Lemma 3. 
Finally, we have: 



Lemma 5: The sequence a at 



lim 



a,N = 1. 




i+ — 

Proof: By Lemma 4, < . Hence, 



N = 1,2, 



is increasing and 



2N 
N 



a N < 4 



The left hand side equals 



2iV 



JV + l 
JV 



2iV 
N 



(j + jfe) 
(1 + 77) 



, from which we get 



2iV 
AT 



jV + l 
N 



< 



2(N+1) 
(N + l) 



which implies that: < o>n+i- To compute limAr^oo a at we observe that 
liniA^oo exp jf [ln(2iV)! — 21n(iV)!] can be computed using the Stirling formula: 



T(x) = 



e x x x 2 



(27T) 



1/2 



1 + 



1 



12a; 



+ 



valid for large x. Thus 



lim oat = 



A/— »oo 



1 r 1 r 

- lim exp — 

4 TV^oo iV 



ln(2iV + l) 2Ar +5 - 2 ln(iV + 



1 

4 iv 



lim e 21n (^) 
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thus completing the proof. ■ 

To find the roots of D' N {9) we solve sin 2 = a^- By Lemma 5 we see that 9jy , s are 
approaching 7r as N increases. Table 2 illustrates this phenomenon numerically. 



Table 2. Roots of D' f 



N 


1 


10 


100 


2 500 


10 000 


root 1 


h 


2.32 


2.80 


3.06 


3.10 


root 2 


h 


3.96 


3.48 


3.27 


3.19 



Let us denote by 6*^-, 9^ the zeros of D' N (0), 9^ € [0, 2n]. Then on the intervals 
In = (@N ~ 2^,^), 1% = (^jv'^jv)' Dn(9) is monotone and thus Dn(9) has an inverse 
which we denote by 0jy and <p 2 N respectively. Now, we observe that the Hilbert space H = 
L 2 (5' 1 ) ~ L 2 ([9^ - 27r, 9jf\) admits the direct sum decomposition: H = £ 2 (/^) © L 2 (I^). 
Moreover, we have: 

L\r N )~L\D N (r N )), n-. m^\{D'NO^ N ){y)\- 1/2 W N {y)), 

where y = D^{9) and i = 1,2. Now we show that r is a unitary isomorphism. 

Lemma 6: Tj is a unitary isomorphism, i = 1,2. 

Proof: 



(V>ilV>2) = / M0)M0) de = 



M<PAv))M<I>n(v)) 



dy 



d9 



dy 



Observe, however, that 



dy 
dO 



{D' N o(f)%)(y)\. Hence, 



(^11^2)^(1*) = / (rMiv) (rMiv) d v = M^rM^iDM 

Jd n (p n ) 



)) 



that ends the proof. ■ 

We will need an explicit form of the inverse of r^. A simple computation yields: 



r- 1 : ^^{D'M^iD^)), 9eP N , <peL 2 (D N (P N )) . 
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Now, we want to define the skew-Hermitian part of the operator p/ 1 ^ on L 2 (Dn). 

It is easier, however, to consider first its push-forward under r: 



1 



— Ti O 

2 1 



2 
1 

-Ti O 

2 

Ti O 



+ 



D' N (9)d9 d9D' N (9) 

1 d d 1 

+ 



D' N (9)d9 d9D' N (9) 
1 (P(D N (9)) + 2 



D' N {8) 

signD' N (9) 



° 7-i V(y) = 

(sign(I^(0)) ^(^)) 1 / 2 

de(j>{D N {6)) 



signD' N (9) 



(si g n(D> N (e))D' N (e)yr< 



(sign(D' N (d))D' N (d)y/* 
d 0( p(D N (9)) = 



5 



Since ^ 



diP Qy y ^ Q~r we need to compute dV Qy V ^ ■ To this end we simply observe that 



D N 0( PN)(y) d( p l N (y) 

dy d^ N (y) 

(<p%°D N )(6) = 9. Hence, 



<W 1 _ 1 
dy ~ D' N (9) ~ (D' N o^ N )(y) ■ 



This implies: 



n o 



d d 
+ 



D' N (9)d9 d9D' N {9) 



° Ti 



(y) = 



d_ 

dy 



yeD N (P N ) . 



We thus have the following important 



Theorem 4: Let X N = \ 



d | d 



D' N (e) de 1 de d' n (6) 



be defined on a dense set S G L 2 {S X ). 



Then X is unitary equivalent to 



L*{D N {I],)) 



dy 



L2(D N (I%)) 



defined on the dense 



space Tl (S n L 2 {P N )) + r 2 (S n L 2 {I 2 N )). 

We can prove a similar statement for Dn(9). We have: 7* o D^{9) o T~ 1 (j)(y) = n o 

lM0)PW(0)l 1/2 <K£;v(0)) = 1^ ° <pKv)\- 1,2 d n o ^(y)l^ ° ^(2/)| 1/2 ^ ° ^)(y) 

= , y e D N (P N ). Hence 

Theorem 5: Let P/v &e i/ie operator of multiplication by the function . Then Pn is 



unitary equivalent to y 



L^D N (li,)) 



®y 



LHD N (I%)) 

Corollary: The pair Q = Dn and P = —X is a conjugate pair on O, where: O = 
T^iSi) ®t~ 1 (S 2 ), Si = {/ G L 2 (D N (P N )),f is absolutely continuous, / = on the 
boundary of -Djv}, i = 1, 2. 
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Remark Observe that X is closed and symmetric on O. We discuss its self-adjoin exten- 
sions below. 

First, however, we determine the spectrum of Dn or to be more precise Td n - By 
the well known theorem on Toeplitz operators [31]: Spec(T Djv ) = essential range of D N = 
range of D N , where the latter is a consequence of continuity. Because X N decomposes 
into a direct sum of operators of the type i-^ on the finite intervals Dn(In) we present 
a brief description of basic features of such operators. For convenience we take a finite 
interval [—1, 1]. The following discussion is based on [32] and all the details can be found 
there. The operator i-^ is known to have a one parameter family of self-adjoint extensions. 
More precisely, let us set i-^ = T a , defined on D(T a ) = {(p : (p is absolutely continuous 
and <f(— 1) = a</?(l)}, where a E C, \a\ = 1. 

To interpret a we solve the eigenvalue problem: = \ip with ip(— 1) = aip(l) = 

e^°V>(l), a = e^°, v?o E [0,2tt). Thus we get: tfj(y) = e lXy and e~ lX = e l(po e iX from 
which it follows that 2A + (po = 2ixn, n E Z. Hence, A = A n = nn — ^p, n E Z. 
The appearance of it is accidental and can easily be removed by rescaling of T a . The 
spectrum of T a is a lattice of physical positions. We would like to point out that there is 
no canonical way of identifying this lattice with the lattice we have started with. It seems 
to be compelling in fact to consider as a real physical space the lattice of eigenvalues of 
T a and the corresponding Hilbert space L 2 (L>at(/^), for example. Thus we would, so to 
speak, avoid the "spectrum doubling " problem by decree. On the other hand this would 
have to be interpreted quantum mechanically as saying that the process of quantization 
amounts in this case to finding an irreducible representation of the Heisenberg commutation 
relations with the additional condition that the spectrum of would be position is discrete. 
This representation would then have a relation to the calculus of difference operators as 
described above, yet, the physical Hilbert space would be only "half" of the Hilbert space 
that naturally carries the action of difference operators. 

To close this section we would like to mention that the optimal discretization has a 
very nice property that limjv^oo Dn(6) = 9. This fact is proven in Appendix B. It is 
only in this case that we get an irreducible representation of the Heisenberg commutation 
relations in the case of the optimal discretization. 

7. SUMMARY AND CONCLUSIONS 

In this paper we have investigated a wide class of discretization schemes for the deriva- 
tive and/or momentum operator defined by (1.3). The optimal subclass of these schemes 
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has been distinguished, that give the best fit to the continuous operators and are sym- 
metric with respect to the space reflections. The particular form of those schemes has 
been determined in Sec. 2. It has been shown that the solution found in [26] is the only 
solution in the form of finite series for the operator X. We also have determined that the 
commutation relation for the annihilation-creation type of operators can be satisfied for a 
large class of difference operators, thus extending the result of [26], if one admits X in a 
form of an infinite series in E. 

Our analysis in Sec. 6 has shown that upon quantization the classically different 
optimal discretization schemes are all unitary equivalent. As a result, the number of points 
used to "delocalize" ^ is unessential. For a fixed number of these points, say 2N + 1, the 
representation of the canonical commutation relations is reducible. Furthermore, by the 
result of Appendix B its irreducible components are unitary equivalent to the case with 
an infinite (N — > oo) number of points. 

Further extensions of this work can be done in several directions. First, one can 
investigate the case of an imaginary shift, as has been mentioned in Sec. 1. Second, more 
general Ansatz for the P and X operators can be analyzed. In particular, we suggest that 
the representation used in [24] can be included into our scheme by considering operator X 
being a nonlinear function of x. Finally, it would be interesting to study simple quantum 
mechanical systems within our approach. 
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Appendices 

Appendix A: Heine symbols 

The Heine symbol defined by (1.7) have pretty interesting behavior. The deformation 
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parameter q > is usually bound to the (complex) region: \q\ < 1 and it is parameterized 
exponentially by a (dimensionless) physical parameter /: 

<Z = exp(^ . (Al) 

The above parametrization allows us to rewrite the Heine symbol (1.7) in the following 
widely used form: 



sm(f ) sinh(^) 
sin(f) 

Function H (P, q) is antisymmetric in P 



sm(i) smh(f) 



H(P,q) = -H(-P,q) . (A3) 
Also, it has the following limits for fixed q: 

lim PYP, q) = , 

P^Q V ' 

hm H(P, q) = 1 , ( A4) 

lim H(P,q) = oo (q^0,l) ■ 
As a function of q the function if (P, q) is invariant under the transformation: 

q^q' = - . (A5) 

For P = 1 it is constant: H(l, q) = 1 and for P < 1 (P > 1) it has maximum (minimum) 
at g = 1 equal to P. The point P = q = 1 is singular and the limits lim^i and limp^i 
do not commute. Also, the limit lim 9 ^o is singular — it is equal for P e [0, 1), equal 1 
for P = 1 and undefined for P > 1. 



Appendix B: A remark about Dn 

Let us consider the function f{9) = 6>, —n < 9 < ir. This function has the Fourier 
expansion: 

oo 

9 = ^ b n sin n6>, 

n=l 

where 

_ 2(-l)" 



n 
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Let us recall that 



n^i n(N + n)\(N - n)\ 

The main idea now is to compare the iV-th partial sum Sn of if (9) with the above sum. 
Clearly, 

JL 2(-l) n+1 
S N (9) = i V -i '- sinnfl. 

' ^ 71 



n=l 



The two sums look alike, the difference is the factor (N-n)\(N+n)} ■ 

We claim that Dn(9) is an approximation to i9, i.e. limjv^oo Dn(9) = i9 in (— tt, it). 
This approximation is a result of a specific method of summation of the Fourier series, 
similar to the Fejer method of arithmetic means. Let us briefly review this subject. We 
refer the reader to [33] for more details. Given an infinite matrix a = [a^- : i,j = 0, 1, . . .] 
and a sequence {S n : n = 0, 1, . . .} we define a new sequence a n , n = 0, 1, ... by 



(TAT = ^ ClNnSn- 



n=0 

We require that a satisfy the following conditions: 

(i) limAr^oo a Nn = 0. 

(ii) £?tv := Yl^Lo l a JVn| exists for every N = 0, 1, . . ., and the set {-Bat : N = 0, 1, . . .} is 
bounded. 

(iii) liniAT^oo Y.n=0 a Nn = 1. 

If the o"at have a limit s then we say that the sequence of partial sums Sn is a — summable 
to the limit s. The most important fact for us is the following: 

Theorem B.l [33]: If a satisfies (i), (ii) and (iii) and if Sn tends to a finite limit s, 
then liniAr^oo ctn = s. 

We claim that Dn are precisely the linear means o~n for an appropriate choice of a. We 
define: 
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■ f l^iSw , if0<n<iV; 

(X Nn — \ (N-n)\(N+n+iy. ' — ' 

0, otherwise. 

It is easy to check that a so defined satisfy (i) and (ii). To see that (iii) is satisfied as well 
we observe that aNn satisfy: 

i (NYf 

£ «»- = (Jy -; ) ,(! y+ „), - s » s w - f^'- 1 * 
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In particular, for n = 1 we obtain: 



A (iV!) 2 N 

1^ aNn ~ (N - l)\(N + 1)! ~ N + 1 ' 

By taking the limit JV -> oo we obtain (iii). To prove our claim we need to show that 
Dn = &n- This follows from (B.l). Finally, applying Theorem B.l we conclude that 
liniTv^oo D N {9) = id. 
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